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1. Introduction 



The relation between classically chaotic systems and the fluctuation properties of 
the corresponding quantum systems has been a focus of intense investigation in recent 
yearsH. We have learned that the quantum counterpart of classically chaotic systems 
have many properties which can be described by random matrix theory (RMT). As a 
consequence, one now typically invokes random matrix ensemble averaging to simplify the 
computation of properties of complex systems. But while the understanding of properties 
of RMT is now quite well established, until recently very little was known about random 
matrix or chaotic Hamiltonians which depend on an external parameter x. The parameter 
dependence can be the strength of an external field, or even the motion of slow variables 
in a complex system treated in the Born-Oppenheimer approximation]^]. The types of 
questions one is interested in are how the properties of such systems are correlated in 
that parameter. In the past few years it has become apparent that a more general class 
of random matrix results exist for parametric Hamiltonians. In a series of papers, Szafer- 
Simons and Altshuler |^ computed the density- density and level velocity correlators in 
a parametric random matrix model, and demonstrated that it described the behavior of 
very different chaotic and disordered systems. By introducing a specific scaling of the 
parameter x, correlation functions were found to fall upon universal curves, independent 
of the underlying properties of the studied systems. Since then, there has been extensive 
studies on universality in other systems, and the class of observables has been extended 
to include wavefunctions and distributions of matrix elements^]. Several further studies 
were based on Ref. including a check that the predictions for wavefunctions hold in 
chaotic systems0, a general theory for scaling^, and the introduction and verification of 
a wide class of universal correlation functions and distributions 0. However, all this work 
has been relegated to one parameter models. The purpose of this paper is to extend the 
previous results to two or more parameters, and to examine the global topological effects 
on universal functions due to the multi-dimensional parameter space. We point out that 
the universality of parametric correlations which emerges on short distance scales can be 
different from 1-parameter model predictions. 



2. Geometric Phase for GOE Matrices 



We consider multiparameter Hamiltonians which exhibit chaos in some region of pa- 
rameter space, such that the fluctuation properties of the quantum Hamiltonian in that 
regime are described by random matrix theory. We will limit the discussion to the gaus- 
sian orthogonal ensemble (GOE), or systems with time-reversal symmetry. As a result, 
the Hamiltonian is a real symmetric matrix, with a non-degenerate eigenvalue spectrum. 
Consider now the general adiabatic variation of the Hamiltonian H['R[t)), where the 
adiabatic eigenstates are defined as 

H{R{t))MR{t)) = E{K{t))M^{t)). (1) 
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For a wavefunction adiabatically transported along some path R, the wavefunction ac- 
quires both a dynamical and a geometric phase, 7|^: 

^(t) = e-(*/^)/o'^*'^"W*'))e*T"W^„(R(t)) (2) 

We are only concerned here with the geometric component of the acquired phase, 7„, 
which contains information about the topology of the parameter space and about the 
presence of diabolical points. For the parallel transport of the adiabatic states ipn around 
a closed circuit C in parameter space, we have 

n, final \ Yn,initial 

) = e'^"^^> (3) 

For real symmetric matrices, the acquired geometric phase for a closed loop is e*'^'-*"^ = 1 
if the path does not enclose a degeneracy (diabolical point), and e*'^'-'"-* = —1 if the path 
does enclose a degeneracy 0. The path can be in any region of parameter space, but if 
the chaotic properties of the Hamiltonian go away for a particular range of parameter, 
one might encounter degeneracies and further complications such as non-abelian gauge 
potentials can develop. We will not consider this here, and hence restrict our discus- 
sion to paths in the classically chaotic regions of the model we study, and equivalently, 
Hamiltonians chosen from the GOE ensemble. 



3. Two Parameter Models 



3.1 Random Matrix Hamiltonians 



Hamiltonian is 10 



We first consider a two parameter realization of the GOE ensemble. A convenient 

(4) 
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(R) = H{x, y) = —= [Hi cos x + H2smx + cos y + H4 sin y] 
v2 

where Ha are independent, N x N, GOE matrices: 





(5) 



-Saf3 i^ik^jl + Sil^jk)- 



The parameter a is related to the average level spacing A via a = Ay/N/TT. H{x,y) is 
clearly GOE at every point {x,y), and has only two relevant moments: 



Hij{x,y) = 



Hij{x,y)Hki{x',y') 



-F{x -x',y- y')i5ikSji + SuSj^ 



(6) 
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with 

F{x,y) = ^ [cos(x) + cos(y)] ^ 1 - -^{x" + y^) + . . . ^ ^7^(0) = i. (7) 

When only one parameter is used Eq. for example keeping y fixed, we have the 

equivalent process: 

a2 

Hij{x)Hki{x') = -i^F{x - x'){SikSji + SiiSjk) , (8) 

where 

Fix) = ^(1 + cos(x)) = cos2(|) ^ 1 - + . . . . (9) 



3.2 A Chaotic Hamiltonian 



We will verify the topological effects we present here in a realistic model. We use 
the Interacting Boson Model (IBM)|]TT|, which is built from scalar (s^) and quadrupole 
= ±2, ±1,0) bosons, which carry angular momentum L = and 2, respectively. 
This model describes the low energy collective excitations of nuclei with even numbers 
of protons and neutrons, the bosons representing paired nucleons. In the consistent-Q 
form||l2|, the Hamiltonian has only two relevant parameters 1] and and in spherical 
tensor notation, is given by: 

H,UV, x) = Eo + vfid + ■ + C3L ■ L, (10) 

Here hd = d) ■ d is the boson number operator, = \/lO[d'^ x d]^^^ is the angular 
momentum, Nb is the boson number defined as half the number of valence nucleons, and 
= djj^s + s^d^ + x[d^ X d]^^^ is the quadrupole operator. Because angular momentum 
is a good quantum number, the parameter C3 adds only an overall constant to Hibm, 
and is unimportant. The typical physical range of the parameters is —^/7/2 < x < 
and < ?7 < 1. The chaotic parameter range of this model has been mapped out in 
detail|]13|, and we consider the quantum properties for this classically chaotic regime. 
The Hamiltonian is diagonalized in the vibrational (or U{5)) basis, with Nj, = 25. The 
universality of one-parameter random matrix theory predictions for wavefunctions and 
distributions has been established for this model and is extensively discussed in Ref. M. 



3.2 Topological Considerations 



The topological properties of the parameter space can be readily taken into account 
using a symplectic 2-form aij and a metric tensor gij. The general Riemannian structure 
of the manifold of quantum states have been defined in Ref. [1^. In the adiabatic basis 
I i/jnCR)) of -?/(R), these are determined by the real and imaginary parts of the quantum 
geometric tensor [Rl p^: 



T 



{ViiJn\{l-\i^n){^n\)\^ji^n) 



9ij + 2^*-'' 
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The antisymmetric tensor aij is related to the element of area of parametric loops we 
study below, and the metric Qij, the measure of distance: 

da = aijdRi A dRj, ds^ = QijdRidRj. (12) 

It is hence possible to incorporate the topological properties of the parameter space into 
a very general formalism. However, as the random matrix model we consider here, in the 
universal regime, is rather simple, we forgo a covariant construction in this article. 



3.4 Universal Scaling for Multiparameter Theories 



Scaling of the parameters is the crucial element in obtaining universality in 1-parameter 
theories. We need to extend this to multi-parameter models. Consider the short distance 
diffusion of the adiabatic energies of a Hamiltonian H{R). From perturbation theory, 
following the arguments of Dyson 0, ^ we have to first order: 



It is convenient to rescale the energies by the mean level spacing £"„ 

1 



2a2 



l-F(R',R)) = (i?:-i?,)(/2; 
-R) DfR) fR'-R). 



r 2N 



(13) 

En/ A, so that : 
(14) 



7r2 ORidR'j 



R'=RJ 



(15) 



Eq. ( [L5[) demonstrates that on short distance scales, the evolution of the adiabatic energies 
of the Hamiltonian resemble a diffusion process, characterized by the diffusion matrix 
Dij(R.). From the above definition, Djj(R) can be related to the autocorrelation of the 
Hamiltonian or the local curvature properties of the energy surfaces: 

2Nd^F{R',R) 



A,(R) 



ViE„(R) ■ V,-E„(R) 



TT 



2 dR'dR', 



(16) 



R'=R 



To obtain universality in correlation functions that depend on the parameters, we must 
rescale the parameters R. In analogy to the 1-dimensional situation, we define: 



R = R 

Here D^^^ is the square root of the diffusion matrix D 



D 



(17) 

J31/2ITJ31/2 N ^ jnodel- 



dependent quantity. Because {SEnY is positive definite, it follows that the diffusion 
matrix Dij(R) must be a positive definite matrix. As a consequence, its square root, 
D-*^/^, is always well defined. Specifically, it is an upper triangular matrix, which can be 
constructed using the Cholesky factorization |jl6|. In terms of R, the energy level diffusion 
is now parameter free: 

Wnf=^' (18) 
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In analogy to the 1- dimensional scaling, we now expect all model dependence to be re- 
moved when we compute correlation functions in terms of R. These are in particular the 
dimension of the Hilbert space N and the short distance behavior of the autocorrelation 
function F{x) = l — cx'^ + ... characterized by the coefficient (s) c. The rescaling by D-*-/^ re- 
moves all this dependence, resulting in parameter free, universal results. We only remark 
here in passing that there is a more general class of processes for which the diffusion is not 
smooth[0, characterized by the short distance behavior of F[x): F{x) = 1 — ex" + ■ ■ ■, 
with < a < 2. In this case, one can proceed as discussed above, but one must define 
the derivatives of F in Eq. (|T^) as fractional derivatives [17, ^. 



In general, there is dependence of the diffusion matrix on its position in parameter 
space, clearly seen in Eq. (p^Gf). A general parametric Hamiltonian H(R) will have an 
autocorrelation -F(R', R) which is not translationally invariant, as is the case for the IBM 
Hamiltonian, so that D is parameter dependent. By construction, our random matrix 
Hamiltonian is translationally invariant, so that F(R', R) = F(R' — R), and as a result 
Dij is independent of R. We find for the random matrix model (Eq. ^): 

Ai(R) = = DSij = —5,j (19) 

In 1-parameter models, D^j D^x = C{0), the scaling introduced in Refs. 0. Note 
that our random matrix model is isotropic, which is certainly not the case for the IBM 
Hamiltonian, which contains bilinear parameter dependence of the form rjx- In that case 
we compute the components of Dij{ri,x) for the energy surfaces En{ri,x) and rescale 
accordingly. 

Consider now the path taken from ^ R. Each of the paths we follow between these 
points can develop a different geometric phase 7„( R). Equivalently, we can consider closed 
circuits which begin and end at a given point, denoted R = 0. Using the two-parameter 
model in Eq. (|), where R = {x,y), we define such a path C as: 

(0, 0) (Xo, 0) ^ (Xo, Yo) (0, Yo) (0, 0) (20) 

Upon traversing this path, the wavefunction will develop a phase according to 

|^n)^e^^"(^)|V^„) = ±|^„). (21) 

From Eq. (|19]), the scaling along the circuit C in the x and y directions is: 

~ [n ^ [n' ^ [n 

Ri = ^—6,,Rj, ^ = \/;^^' y^i^y- ^^^^ 

The relevant quantity is the area of the loop, A = AxAy = XqYq. This quantity can be 
made more useful for comparisons to parametric areas in other Hamiltonians by defining 
the scaled area: 



A = AxAy = ^DxxDyyXoYo. (23) 

As a check of our results, we plot in Fig. 1 the ratio of the computed value of A to the 
predicted result in Eqs. (p2D , (p3D , using the Hamiltonian of Eq. (4), for varying sizes of 
the parametric loop. We take 100 points around the loop C, and compute Dij{x,y) by 
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Figure 1: Comparison of computed scaled area, A — ^J^^^D^XqYq, of the parametric 
loop C, for the Hamiltonian of Eq. (4), to the theoretical result from Eqs. (22)-(23). The 
dashed line is the expected result. 
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averaging over the middle third of the energy surfaces En{x,y). The agreement is quite 
good. (There is a shght systematic shift of the results which seems to be due to the 
method used to unfold the energy spectrum.) 

The regime of universality is roughly defined by the scale M 



R 



< 



1. (24) 



Analogously, the regime of universal behavior for topological phase accumulation is 

i ~ 1. (25) 

We will see in the next section, that the area of the loop does not have to be large 
before one finds significant effects, and saturation of these effects already occur within 
the universal regime. 

Our previous work in Refs. |^, ^] examined the universal statistical decorrelations 
of wavefunction related observables in a one parameter formalism. In Fig. 2, a typical 
universal result is shown for the distribution of wavefunction overlaps, P(m;R), with 
u = ('?/'n(R') I ^n(O)), for a 1-parameter path. (The vertical axis in all plots of P(m; R) 
is rescaled to place the maximum near unity.) The solid histogram is the result from 
random matrix theory, and the dashed histogram is the result from the IBM, in particular 
J'^ = lO"*" states 0. Using Nfy = 25, this corresponds to a dimension of 211 states, of 
which only the middle third of the eigenstates are used. The distribution begins as a 
delta function at R = 0, and develops into a Porter-Thomas distribution for R ^ 1. The 
Porter-Thomas limit is non-universal in the sense that the distribution depends explicitly 
on the dimension of the space. A one-dimensional path in the 2-dimensional parameter 
space (Eq. 9) has the same universal behavior when plotted as a function of the scaled 
distance R = Da ■ R. However, if we now allow paths in two or more dimensions, the 
shape of this distribution depends on how one reaches the position R. 

Consider the same distributions when we traverse the path of Eq. (pO|) of area A. 
When we transport A^ eigenstates around the loop, a certain number of them will develop 
a phase 7n(C) = ti", which we denote A/'(7r). The remaining A/'(0) states do not acquire a 
phase, where A^ = A/'(0) +^/{tt). The fraction of states which acquire a Berry's phase is 
defined as: 

^ ~ J\f{0)+M{7ry ^^^^ 

As mentioned above, because a non-zero value for //{it) indicates that paths followed 
by certain eigenstates enclose a diabolical point, this measure is clearly related to the 
number of diabolical points enclosed by the path of area A. In Fig. 3, the wavefunction 
overlap distribution function, P{u; R), is shown at various points along the square circuit 
C. The solid histogram corresponds to the random matrix predictions, and the dots 
to the results from the IBM. Initially, the distribution is a delta function centered at 
{x, y) = (0, 0) . As the separation increases, the decorrelation is similar to the results of 
Fig. 2. But as the circuit returns to the origin, the distribution bifurcates into two distinct 
distributions, corresponding to the existence of the Berry's phase of ±1. As the circuit is 
closed, the distribution does not go back to its original form, but is now described by two 
delta functions located at ±1 of equal magnitudes (within numerical fluctuations). The 
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Figure 2: One-parameter universality of the distribution of matrix elements P(m;R), 
where u = {ipn(R^) \ "^nlO)). The solid histogram is random matrix prediction while the 
dotted histogram is the observed result from the Interacting Boson Model 0. 
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fraction / is at its saturation value of 50% for this loop, as seen in Fig. 3(f). Hence, the 
short distance universality that occurs in multiparameter systems is not simply a function 
of the metric distance between the two points and R. Given only the distance R, the 
distribution P{u; R) is not well defined, since one might be close to either Fig. 3(a) or 3(f), 
which represent the same point. However, both are universal results when one follows a 
path of the same scaled area. In this case, the universal function P{u) must also identify 
the area subtended by the path. Figure 4 shows a similar collection of distributions in the 
case of a smaller parametric square. As can be seen here, only a small fraction of states 
subtend a diabolical point. It is clear that quantities which are sensitive to the phase of 
the wavefunctions will generally be modified in multi-parameter theories. 



4. Universal Density of Diabolical Points 



The effect we have described in the previous section is due to diabolical points which 
are enclosed by the parametric paths. We now consider the density of these diabolical 
points. It has been shown that a necessary condition for the occurrence of a Berry's 
phase of 7(C) = vr for a system transported around a closed loop C is the existence of 
a diabolical point within the loop]^, |T^. Hence the fraction / of states in Fig. 3(f) and 
4(d) that have 7(C) = tt should be a measure of the number of diabolical points enclosed 
in the area A, or equivalently, the density of diabolical points. We note that we assume 
that as the area increases, the occurrence of a new diabolical point is responsible for the 
phase change of only one eigenfunction. For loops of small area, which is the situation in 
the universal regime, this is a reasonable assumption. In Fig. 5 (top) we plot the fraction 
of states / which enclose a diabolical point as a function of the scaled area. There are 
three sets of results in the figure. The crosses are the result of varying the size of the 
loop, Xq = Yq E (0,0.25], with = 200, so that the scaled area varies from to 1.27. 
The boxes correspond to varying the dimension of the matrices from = 50 — 300, for 
fixed area with Xq = 10 = 0.18. Finally, the open circles correspond to parametric loops 
in the Interacting Boson Model, at fixed dimension. As is clearly seen the behavior is 
statistically equivalent for all the results, demonstrating universality of the results. The 
general features of these results are the following. We see that the fraction / increases 
linearly with scaled area, and that statistical saturation of 50% occurs for A ~ 1. (In 
the Hamiltonian (4), the relation to the area is A = NA/ir"^ (see Fig. 1) while in the 
IBM, A oc NA.) We expect saturation when there are many diabolical points, so that 
statistically the probability of having a phase change for an arbitrary wavefunction is 1/2. 
If we define the number of diabolical points enclosed by the circuit C as n{C) = Nf{C), 
we conclude that in the non-saturation regime, 

f{C) = cNA, n{C) = cN^A, (27) 

where c is determined from the slope of Fig. 5 (bottom). In the saturation regime, 

f{C) = I, n{C) = ^. (28) 
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Figure 3: Modification of 1-parameter universality for a 2 parameter square of side Xo = 
Yq = 0.32, with N=200. The distribution P{u; R) is shown at the points (X,Y) = (a) 
(0,0), (b) (0.16,0), (c) (0.32,0.32), (d) (0,0.32), (e) (0,0.16), (f) (0,0) (after traversal 
of loop). Both (f) and (a) represent the same point, showing that universality is path 
dependent. The dotted histogram corresponds to results from the IBM, on a parametric 
loop of similar scaled area. 
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Figure 4: Same as Fig. 3, but for a loop of smaller area. In this case, Xq — Yq — 0.03, and 
N = 200. The figures correspond to {X, F) = ( in units of Xq) (a) (0, 0), (b) (1.0, 1.0), (c) 
(0.8, 1.0), (d) (0, 0) (after traversal of loop). The smaller loop area contains fewer diabolical 
points and hence has fewer states which split in (d). The dashed histogram corresponds 
to the results from the IBM, which are plotted slightly offset for better contrast. 
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When the area is large, / is no longer a good measure of the number of diabolical points, 
so the dependence of n{C) and f{C) with N is no longer meaningful. In Fig. 5 (bottom), 
we plot the small area behavior of the fraction /. The statistical errors are computed 
by computing the fraction / and the scaled area A for a number of realizations of the 
parametric Hamiltonians. The results for the chaotic system (IBM), the open circles, are 
not averaged over, but come from a single parametric loops. The general result is that 
the number of diabolical points grows quadratically with the dimension of the matrix. 
(Saturation effects occur for larger areas). From the results of Fig. 5, we conclude that 
the fractional density of diabolical points, p, in terms of scaled area, is a universal constant 
for chaotic systems, given by 

p = = 0.94 ±0.10 (29) 

The quantity similar to n{C) was examined recently in a very interesting study of avoided 



level crossings [T^, using the same two-parameter GOE model. In that study, a similar 
scaling in was found for the number of diabolical points, although their results were 
not universal, and do not completely agree with the results of this study. The overall 
scaling agreement is encouraging, and we believe the discrepancy is due to the different 
methods used to count diabolical points. 



5. Conclusions 



We have developed a formalism to compute universal parametric correlations in multi- 
parameter systems, explicitly demonstrating universality in chaotic systems. The scaling 
in achieved through the square root of the diffusion matrix Dij. We have also found that 
universal results for two-point correlation functions and distributions, which are thought 
to be only a function of the separation of the points in parameter space, can be path 
dependent in more than one parameter dimension. Indeed, one can obtain distinct results 
for the universal correlations by reaching the same final point by different paths. So short 
distance correlations are not a simple function of the metric distance between the points 
in parameter space, but are also a function of the shape of the path taken between the 
points. This is true for all quantities which are sensitive to the phase of the wavefunctions. 
Whether one can actually measure interference effects in a manner analogous to measures 
of Berry's phase is an interesting question to explore. Quantities which are not sensitive 
to the phase, will not suffer from this ambiguity. Finally, we have verified that the density 
of diabolical points grows quadratically with the dimension of the matrix, and used that 
to determine a universal measure of the density of diabolical points in chaotic systems. 

This work was supported in part by the Department of Energy Grant DE-FG02- 
91ER40608. 
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Figure 5: Universality of the fractional number of diabolical points /. (top) Fraction 
of wavefunctions enclosing a diabolical point as a function of the scaled area of the loop 
C, calculated at (i) fixed dimension = 200, with the size of the loop varying over the 
range < Xq = Yq < 0.25 (crosses); (ii) fixed loop area with Xq = Yq = 0.18, and the 
dimension varied from N — 25 — 300 (boxes); (iii) fixed dimension in the chaotic region of 
the IBM parameter space, for several loop sizes (circles), (bottom) Linear behavior of the 
small area behavior of /. From this we find the the density of diabolical points increases 
quadratically with the dimension of the matrix. 
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